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ABSTRACT 

We compare the accuracy of published formulae that transform the linear perturbation 
theory power spectrum into the nonlinear regime against the results of an ensemble of 
large N-body simulations, paying attention to scales on which the density fluctuations 
are linear and quasi-linear. The inverse transformation to obtain the linear power 
spectrum is applied to the APM Galaxy Survey power spectrum measured by Baugh 
& Efstathiou (1993). The resulting linear spectrum is used to generate the initial 
density fluctuations in an N-body simulation, which is evolved to match the measured 
APM amplitude on large scales. We find very good agreement between the final power 
spectrum of the simulation and the measured APM power spectrum. However, the 
higher moments for the particle distribution only match the ones recovered from the 
APM Survey on large scales, R > 10ft.~^Mpc. On small scales, R < 10ft.~^Mpc, the 
APM estimations give smaller amplitudes, indicating that non-gravitational effects, 
such as biasing, are important on those scales. Our approach can be used to constrain 
a model of how light from galaxies traces the underlying mass distribution. 
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1 INTRODUCTION 

The growth of density fluctuations can be followed accu- 
rately using linear perturbation theory only when the den- 
sity contrast on a given scale is much smaller than unity 
Sp/p <^ 1 (see for example Peebles 1980). As the fluctua- 
tions enter the mildly nonlinear regime, Sp/p ~ 1, analytic 
approximations or numerical simulations have to be used to 
follow the evolution of the density field, except in cases with 
idealised geometry {e.g. Bertschinger 1985). 

Several comparisons between the results of N-body sim- 
ulations and the predictions of higher order perturbation 
theory have been made recently {e.g. Buchert, Melott & 
Weiss, 1993, Jain & Bertschinger 1994). Baugh & Efstathiou 
(1994b) demonstrated that Eulerian second order perturba- 
tion theory gives a good approximation to the evolution of 
the power spectrum in Standard Cold Dark Matter (flat 
universe, £7 = 1, with h — 0.5^: hereafter SCDM) down to 
scales for which the variance ^2 ~ 1 correctly predicting a 
transfer of power from large to small scales. 

An alternative approach for following the evolution of 
the density field into the nonlinear regime has been adopted 



The Hubble constant is given by Ho = lOO/ikms-^Mpc" 



by Hamilton et aL(1991), who deduced a transformation be- 
tween the linear volume averaged correlation function and 
the nonlinear correlation function, with the functional form 
calibrated against the results of numerical simulations. Pea- 
cock and Dodds (1994 - hereafter PD) proposed that a simi- 
lar form of the transformation could be applied to the power 
spectrum of density fluctuations. Jain, Mo and White (1995 
- JMW) suggested a correction that improves the perfor- 
mance of the Hamilton et al. and PD formulae when the 
effective slope of the power spectrum, P{k) cx fc" has the 
value n < ~1. 

In Section 2 of this Letter, we compare the predictions 
of the formulae of PD and JMW for the evolution of the 
shape of the power spectrum in an SCDM universe against 
the results of an ensemble of large N-body simulations. 

We use these formulae to compute the linear power 
spectrum corresponding to the power spectrum measured for 
APM Survey galaxies by Baugh & Efstathiou (1993, 1994a) 
in Section 3. We use this linear power spectrum to generate 
the initial conditions in an N-body simulation. The clus- 
tering in the evolved particle distribution is then compared 
with the measurements of the APM power spectrum and the 
higher order moments of counts in cells (Gaztafiaga 1994). 
This approach provides a test of our ideas about structure 
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Figure 1. The ratio of the power spectrum at epoch to the 
power spectrum at some reference epoch ao = 0.4 (where a = 1 
corresponds to z = 0), with the linear growth scaled out. The 
points show the results from an ensemble of large N-body sim- 
ulations. The lines show the predictions of the PD and JMW 
formulae for the linear to nonlinear transformation. The arrow in 
the middle panel shows the size of the simulation box used in the 
largest volume run of Jain et al.. In the bottom panel, the lines 
show the fit quoted in the text. 



formation, such as the Gaussianity of the initial conditions 
and the biasing of the galaxy distribution relative to the 
mass distribution. 



2 EVOLUTION OF THE CDM POWER 
SPECTRUM 

In this Section we shall examine the CDM power spectrum 
as an example of a scale dependent spectrum. The nonlinear 
growth of density fluctuations can be studied by comparing 
the shape of the power spectrum of the fluctuations at dif- 
ferent epochs. Following Baugh & Efstathiou (1994b), we 
take the ratio of the power spectrum at expansion factor 
corresponding to time ti, to the power spectrum at some 
reference epoch, set by ao and time to, scaling out the linear 
growth factor oc a?: 



R{k,U;to) 



P{k,U) 
P(k,to)' 



(1) 



This approach takes into account any inaccuracies in the ini- 
tial realisation of the theoretical power spectrum (see Baugh 
& Efstathiou 1994b). We have found that an expansion fac- 
tor of about 3 is necessary to remove transients arising from 



the Zeldovich approximation that is used to set up the ini- 
tial fluctuations (see Baugh, Gaztafiaga & Efstathiou 1995, 
hereafter BGE95), hence our choice of a = 0.4 as a reference 
epoch. 

The general form of the transformation of linear to non- 
linear power proposed by PD is given by 



kL = [l + A%L{k^ 



-1/3, 



(2) 
(3) 



where the subscripts L and NL refer to linear and non lin- 
ear respectively and A(fc) = 4tt k^ P{k)/ {2tv)^ is the frac- 
tional variance of the density field in bins of In k. The JMW 
correction requires making a shift in this universal relation: 



A\kNL)/b{n) = /jvL[Ai(fci)/6(n)], 



(4) 



where the factor b{n) is a function of the spectral index of 
the power spectrum. JMW obtained the form of this fac- 
tor b{n) by requiring that their transformation reproduced 
the results of simulations with scale free initial conditions. 
They then make the assumption that these results can be 
applied to the case of fiuctuations that do not have a scale 
free initial power spectrum, such as CDM. In this instance, 
an effective spectral index n^, is defined as the local slope 
of the power spectrum on the scale at which the variance in 
density fluctuations is unity. 

Using the flts given by equation (23) of PD and equa- 
tions 5(b) and 7(a) of JMW, we can generate a set of evolved 
power spectra at different epochs. We set a = 1 at ag = 1 
and choose a — 0.40 to be the reference epoch. We interpo- 
late over the power spectra generated at a = 0.50, 0.59, 1.00 
to form the ratio in equation (^). The ratios predicted by 
the formulae are shown as the solid lines in Figure |^. Also 
plotted are the same ratios calculated from an ensemble of 5 
N-body simulations, with a box size of 378/1" ^ Mpc and 126^ 
particles. The simulations were run with the P^M code de- 
scribed by Efstathiou et aZ.(1985). The errors are the disper- 
sion in the ratios over the five simulations in the ensemble. 
We have made no correction to the N-body results to take 
into account aliasing arising from assignment of particles to 
the FFT grid, as these corrections largely cancel out when 
the ratio of power spectra is taken, due to the relatively small 
change in spectral index between different epochs. The cor- 
rection for particle discreteness is uncertain (BGE95) and is 
small due to the large number of particles used. These ef- 
fects make some contribution to shape of the N-body curves 
at very high wavenumbers around the particle Nyquist fre- 
quency, k > 1/iMpc^^, but do not affect the comparison on 
large scales. 

In Figure ^, we show a fit for the evolution of the power 
spectrum which gives a good match to the results of the large 
box N-body simulations used in this paper. We have used the 
same type of fit as JMW, with b(ne) = 1.16[(3 -I- ne)/3]° '^ 
where rie is the effective spectral index at each epoch as 
defined by JMW, with /nl given by 



/jVL (x) = X 



1 + ax + bx^ + cx"^ + da;^'^ + ex' 



1 + [e/(11.68)2] 



1/2 



(5) 



We find a = 0.598, b = -2.390, c = 8.360, d = -9.010, 
e = 2.895 and / = —0.424 by matching the power spectrum 
in the simulations at a = 1, i.e. erg — 1. The accuracy of the 
fit is better than 5% over the range 0.02 < fc < 1.0/iMpc"^ 
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Figure 2. Open circles sliow tiie APM power spectrum with Icr 
errors. The lines show the linear spectrum predicted by the JMW 
(solid) and PD (dotted) formulae. Symbols show P{k) measured 
from a N-body simulation with the initial conditions set by the 
JMW (filled circles) and PD (open squares) linear spectra. 

Note that as our simulations do not have the resolution to 
probe the highly nonlinear regime, where we have forced the 
fit to have the asymptotic form fNhix) = ll.QSx'^^^ when 
X —> oo, as used by PD. Thus our fit does not necessarily 
perform well on small scales, where JMW results are more 
reliable. 

Whilst the agreement of the JMW formula with the 
simulation results on large scales is within the quoted 20% 
accuracy, the transformation performs less well than either 
PD or the fit given in this paper. This is mainly the result of 
an overprediction of the power at early epochs. This could 
be due to a number of reasons. JMW used high resolution 
simulations in order to examine the behaviour of the power 
spectrum in the highly nonlinear regime. This was achieved 
by using a relatively small box size, as indicated by the arrow 
in the middle panel of Figure |l] Hence, Fourier modes of the 
density field around k ~ 0.2/iMpc~^ do not have modes on 
larger scales to couple to, with the result that the nonlinear 
evolution on these scales cannot be followed accurately. 



3 N-BODY REALISATIONS OF CLUSTERING 
IN THE APM SURVEY 

The PD and JMW formulae for the inverse transformation 
from nonlinear to linear power can be used to obtain the 
linear power spectrum that corresponds to the APM Sur- 
vey power spectrum (APMPK) measured by Baugh & Efs- 
tathiou (1993, 1994a). We make the assumption that there 
is no bias between the galaxy distribution and the mass, i.e. 
that the light traces the mass. 

Figure ^ shows the APMPK as open circles with Icr er- 
rors obtained by averaging over the APM Survey split up 




k 



Figure 3. The logarithmic slope of P{k) as a function of 
wavenumber estimated from the measured APMPK (circles with 
errors) and from the N-body simulation started from initial con- 
ditions predicted by the JMW formula (filled squares). The er- 
rorbars are computed using the Icr scatter on the power spec- 
trum. The long (short) dashed curve shows the slope in the linear 
r = nh = 0.5 (r = 0.2) CDM model. The solid (dotted) line 
shows the fit to the measured (linear) APMPK. 



into four zones. The solid line shows the linear power spec- 
trum predicted by JMW, which has a slope of n ~ —2, 
and the dotted line shows that obtained with the PD for- 
mula. The kink in the linear spectra at fe > l/iMpc"'^ is 
an artefact due to the form of the transformation formulae. 
The simulations that we use are designed to a have parti- 
cle Nyquist frequency at shorter wavenumbers than this and 
are not affected by this feature. The two linear spectra are 
quite different. We find that the simulations evolved from 
the JMW formula give better agreement with the shape of 
the APMPK, confirming that this transformation works best 
when the power spectrum is nearly scale free. The filled cir- 
cles show the power spectrum of a P^M simulation with 
160"^ particles and a box size of 440/i"^Mpc. The simula- 
tion has been evolved over ~ 6 expansion factors from the 
initial conditions to match the amplitude of the variance 
in spheres of radius 30/i~^Mpc given by the APMPK. Also 
shown is the result of a smaller simulation, using the PD lin- 
ear power spectrum (open squares), which shows an excess 
of small scale power compared with the APMPK. PD note 
that their formula is not expected to work well for such a 
steep spectrum (see also Figure 2 of JMW). We note that 
the linear spectrum predicted by our fit is closer to that of 
PD rather than that of JMW; hence a transformation of 
the type given in equation ^ which is calibrated against the 
results of CDM simulations does not work well when used 
with an almost scale free spectrum. 

Figure 3 shows n = din P(fe)/dln fe, the logarthmic 
slope, estimated from the measured APMPK as a function of 
wavenumber (open circles). The dotted line shows the slope 



© 0000 RAS, MNRAS 000, 000-000 



4 CM. Baugh and E. Gaztanaga 



of a fit (see below) to the linear APMPK using the JMW 
formula. The filled squares show the slope of the power spec- 
trum after evolution in the N-body simulation, which repro- 
duces closely the form of the observed APMPK. This Figure 
clearly shows that the feature in the observed APMPK at 
k ~ 0.3Mpc~^ is produced by nonlinear evolution, as it is 
not present in the linear APMPK. Nonlinear evolution in 
CDM-like models produces similar behaviour, but is insuf- 
ficient to explain fully the shape of the measured APMPK 
(Baugh & Efstathiou 1994b); the linear APM spectrum is 
too steep to be fitted by any CDM-like model, expressed as 
a function of the parameter V — flh (Efstathiou et a/. 1992), 
which are shown as dashed lines in Figure ^. Efstathiou et 
al. 1990 have shown that it is difficult to match both the 
large and small scale shape of the angular correlation func- 
tion measured in the APM Survey with CDM-like models. 
We note however, that the shape of the linear APM power 
spectrum is similar to that predicted in Mixed Dark Matter 
models (Efstathiou 1995). 

To produce the perturbation theory predictions in Sec- 
tion 4, we fit the estimated P{k) to a simple parametric 
shape: 



P(fe) 



C fc" 



[1 + (fc/fcc)2]' 



(6) 



with a minimum fit using the estimated errors. Given the 
additional uncertainties from the APM selection function, 
evolution of clustering and value of il (see Gaztanaga 1995) 
we use 2-sigma errors in the APMPK as the error estimation 
in this fit. The results are shown in Figure ^. For the non- 
linear APMPK a fit to the whole range of k gives (solid 
line): C ~ 9.5 x 10^ kc ~ 0.03/i"^Mpc, a ~ 1 and & ~ 1.2. 
The fit to the linear JMW P(k) (dotted line in Figure ||) 
is restricted to < 0.6 and gives: C ~ 7.0 X 10^ kc ~ 
0.05/i~^Mpc, a ~ 1 and 6 ~ 1.6. The reduced is much 
lower for the fit to the hnear curve than for the nonlinear 
P{k), showing the difficulty in reproducing the interesting 
features in the APMPK around k > 0.3Mpc~^ with a simple 
parametric form. It is interesting to note that although a is 
a free parameter the best fit gives in both cases a ~ 1, as 
predicted by the Harrison-Zeldovich spectrum. 



4 HIGHER ORDER MOMENTS 

We next evaluate the higher order moments in the N-body 
simulation with the same power spectrum as the APM 
survey. We use the counts in spherical cells of radius R, 
to estimate the volume averaged J-order correlation func- 
tions S,j{R), as described in BGE95. We concentrate on the 
higher order moments in terms of the hierarchical ampli- 
tudes Sj = 0/^2 ■ These quantities can be predicted in 
perturbation theory for models with Gaussian initial fluctu- 
ations which evolve only under gravity. Bernardeau (1994) 
has estimated Sj for the case of a spherical (top-hat win- 
dow) cell, which are given in terms of J-order logarithmic 
derivaties fj ol (^^{R) in the initial conditions. These predic- 
tions have been tested up to J = 10 in N-body simulations 
(Gaztafiaga & Baugh 1995, BGE95), showing a very good 
agreement for scales where ^2 i$ 1- 

In Figure ^ we compare the perturbation theory predic- 
tions for Sj with J = 3, 4, 5 using both the linear (solid line) 




4 6 8 10 
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Figure 4. The open circles show the APM estimations for S3, 
S4 and 55. The closed triangles show the corresponding moments 
measured in a large N-body simulation started with the linear 
APM P(k). The solid (dashed) line shows the perturbation theory 
prediction for Sj using linear (non-linear) APM power spectrum. 



and non-linear (dashed line) shape of ^2 estimated from the 
P(k) fits above. Although we can see in Figure ^ that these 
predictions are quite different at small scales, they agree 
well on scales where C2 ~ 1- On comparing with the evolved 
results from the APM like N-body simulations (closed trian- 
gles in Figure ^) , we find a good agreement with perturba- 
tion theory predictions for scales where ^2 J; 1) expected. 



5 DISCUSSION 

The modified transformation from the linear to nonlinear 
regime introduced by JMW seems to work well for the case of 
the linear APM spectrum presented in Section 3. However, 
the transformation is less successful in predicting the de- 
tailed nonlinear behaviour for the SCDM model. The trans- 
fer of power in the midly nonlinear regime is of the same 
order of magnitude as the quoted accuracy of the JMW fit. 
Whilst we have shown that an improved fit to the behaviour 
of CDM on large scales is possible by constraining the for- 
mula to match the results of simulations with a large vol- 
ume, there appears to be a fundamental problem in getting 
a truly universal fit that works equally well for scale free 
and scale dependent spectra, with comparable accuracy on 
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large scales to that possible with second order perturbation 
theory (Baugh & Efstathiou 1994b). 

In previous comparisons of perturbation theory predic- 
tions for Sj{R) with clustering in the APM Survey (e.g. 
Gaztaiiaga 1994, Gaztanaga & Frieman 1994, BGE95) the 
non-linear shape of P{k) or used to calculate the 

perturbation theory predictions. This assumes that on large 
scales the observed shape is not affected by non-linear bias- 
ing or by non-linear gravitational effects. Here we drop the 
latter assumption and use the linear P{k), as required in 
perturbation theory. The comparison of these improved pre- 
dictions with the estimations from the APM Survey (open 
symbols) is shown in Figure ^. 

The APM estimations are the ones presented by 
Gaztaiiaga (1995) which result from using a simple scaling 
law to model the projection effects. Although there are some 
potential problems with this modelling (Bernardeau 1995), 
we believe that these results are accurate (see Gaztaiiaga & 
Baugh 1996). The APM amplitudes agree quite well with 
the improved predictions on scales R > 3/i~^Mpc. This is 
surprising as one would rather expect to find an agreement 
with the fully evolved N-body results, which do not match 
these predictions at small scales. Thus, our analysis indi- 
cates the failure of at least one of the hypotheses we have 
used. 

We have used fl = 1. This is not very important for the 
2D to 3D inversion of the APM correlations, as changing 
the cosmology only alters the overall clustering amplitude 
slightly and not the shape of the correlations (see Baugh & 
Efstathiou 1993, Gaztanaga 1995). In the N-body models, 
a different value of Q would change both the infered initial 
P{k) and the final Sj in Figure W. We have run some test 
models and find that for the APMPK the fitting formulas 
do not seem to work that well for f2 < 1. We find that the 
spectral index is predicted to be slightly more negative on 
small scales when fl < 1, compared with the linear power 
spectrum obtained for Q — 1. This means slightly larger 
perturbation theory predictions for Sj and also more non- 
linear evolution at high k. We find nevertheless very little 
difference for the final values of Sj for different values of Q. 

In our APM-like simulations we have assumed Gaussian 
initial conditions and in order to infer the linear mass power 
spectrum from the measured galaxy power spectrum we have 
also assumed that there is no bias between the fluctuations 
in the galaxies and in the underlying density field. The hier- 
archical scaling of the higher order moments measured from 
the APM Survey (Gaztaiiaga 1994) suggests that there is 
no relative bias between mass and light on large scales. In 
addition it is unlikely that the scaling could be produced by 
non-gaussianities or a particular biasing prescription that 
happens to mimic gravitational growth (Gaztafiaga & Frie- 
man 1994). However, the disagreement shown in Figure ^ 
indicates that either of these assumptions fails on scales 
R < 10ft"^Mpc. Some form of non-gravitational effect or 
small scale dependent biasing is necessary. Non-gaussian ini- 
tial conditions with Sj < on small scales, could account 
for the smaller values of Sj in the measurements. Biasing 
would also alter the galaxy amplitudes Sj both directly and 
through the change of the underlying mass power spectrum 
which will lead to a different prediction for the linear mass 
power spectrum. 
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